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ON THE EFFECTS OF MONOENERGETIC ELECTRONS IN A 


STEADY-STATE, LOW PRESSURE PLASMA 
by Nelson L. Milder 
Lewis Research Center 

SUMMARY 

This paper describes a steady-state free-fall theory of a one-dimensional, low- 
pressure, beam- generated plasma in a plane diode configuration. The motivation for 
this work derives from the fact that the low pressure plasma discharge encountered in 
electron-bombardment thrusters is believed to develop an electron density composed of 
both monoenergetic and Maxwellian electrons. It is of interest to determine the effects 
of such a distribution on plasma potential, plasma-sheath boundaries, ion current den- 
sity, and mean ion energy. For this reason the earlier theory of the positive column 
introduced by Tonks and Langmuir was expanded to include an electron density consist- 
ing of a superposition of cathode emitted, collisionally scattered monoenergetic (pri- 
mary) electrons upon a Maxwellian distribution. Cases are considered in which the pri- 
mary electron space charge contributes to plasma neutralization as well as ion produc- 
tion. In such cases, the potential distribution, ion current density, and mean ion energy 
in the plasma differ from the primary-free case. 

The presence of a monoenergetic electron space charge density made the usual 
Tonks- Langmuir plasma-sheath boundary condition inapplicable. (The Tonks- Langmuir 
condition is that the reciprocal of the spatial variation of potential is zero at the sheath. ) 
Instead, the Bohm sheath criterion was used to replace the Tonks- Langmuir boundary 
condition by one that depends on primary electron collision phenomena. Plasma poten- 
tial depended on the primary electron density, ionization frequency, and to a lesser ex- 
tent, beam scattering phenomena. 

Approximate models of the plasma discharge are considered in which the primary 
electrons contribute to ionization but not to the space charge density. With the appro- 
priate choice of parameters, potential profiles obtained from the approximate models 
coincide with profiles from a more complete model over a wide range. The solutions 
based on the approximate models are more easily obtained and can be used to investi- 
gate appropriate regions of low pressure plasmas. 


INTRODUCTION 


This paper describes a steady-state theory of a beam- generated plasma in which 
beam electron scattering by collisions is included. The motivation for this work derives 
from the fact that the low-pressure plasma discharge encountered in electron- 
bombardment thrusters is believed to develop an electron density composed of both 
monoenergetic primary electrons and Maxwellian electrons (ref. 1). It is of interest to 
determine the effects of such a distribution on plasma phenomena. 

The analysis considers a one -dimensional configuration associated with a plane 
parallel diode system. The diode consists of a plasma discharge region between a 
cathode emitter and an anode. In previous work on beam- generated plasma in a plane 
symmetric configuration (ref. 2), the electron beam flowed parallel to and between two 
grounded plane electrodes. These electrons were presumed to suffer no scattering col- 
lisions and were the principle source of ionization. The present study considers effects 
of an electron distribution consisting of cathode -emitted monoenergetic primary elec- 
trons superimposed upon a Maxwellian electron distribution. The primaries are allowed 
to suffer nonionizing as well as ionizing collisions. The details of such collisions need 
not be known explicitly. 

For the case in which primaries contribute to the space charge density in the 
plasma, a low pressure approximation is used. In essence, this means that the primary 
electron mean free path X is of the order of or greater than system dimensions. The 
primary electrons in the plasma are assumed to follow an exponential spatial distribu- 
tion ~e -Z// ^ in their flow from cathode to anode. Ion generation is taken to be propor- 
tional to contributions from both primary and Maxwellian electrons. 

The theory is based upon an extension of the one-dimensional free-fall theory of 
Tonks and Langmuir (ref. 3). Starting with the one -dimensional Poisson equation, an 
integral equation is derived in reference 3, which is presumed to be valid in both the 
plasma and collisionless sheath regions of the discharge. This equation is the Tonks- 
Langmuir plasma-sheath equation. Solutions to this equation yield the plasma potential 
distribution. The original theory, in which second order terms could be neglected, re- 
sulted in solutions for the potential and current distributions in plane, cylindrical, and 
spherically symmetric low-pressure plasma arcs. Interest in this theory was increased 
by the work of Harrison and Thompson (ref. 4) on a low pressure, plane symmetric dis- 
charge. Their paper stimulated further work by several authors (refs. 5 to 9). Recent 
efforts extend the work of Dunn and Self (ref. 2) to a cylindrically symmetric beam- 
generated plasma (ref. 10) and provide experimental verification of the collisionless 
sheath theory (ref. 11). 

The present calculations yield information on the effects of different modes of ion 
generation on the plasma potential distribution and ion current density in the plasma due 
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to the presence of monoenergetic electrons. In addition, the mean ion energy and Bohm 
sheath criterion are investigated in the context of the theory. The presentation begins 
with a description of the physical system to be considered The analysis incorporates 
an electron distribution consisting of primary beam electrons plus Maxwellian electrons. 
In the limit of vanishingly small ratio of Debye length to characteristic plasma length, 
this equation reduces to the so-called plasma equation in which charge neutrality holds. 
Solutions to this plasma equation are investigated. 


PHYSICAL SYSTEM 

A schematic diagram of the diode configuration to be studied is given in figure 1. 

In the figure the plasma potential at potential maximum is assumed zero. Electrons 
emitted from the cathode gain kinetic energy in the cathode sheath and arrive at the 
potential maximum with an energy V , where V c is the cathode potential. The dashed 
line represents the potential variation in the sheath regions. As a result of ion produc- 
tion and randomizing processes, the two kinds of electrons presumed to exist in the 
plasma discharge region are cathode emitted monoenergetic electrons (primaries) and 
Maxwellian electrons at uniform temperature T. Primaries suffer collisions (ioniza- 
tion, excitation, scattering, etc. ) when they traverse the discharge. Once a primary 
electron suffers an energy absorbing or scattering collision, it is assumed to be re- 
moved from the beam, so that the primaries follow an exponentially decaying spatial 
distribution. 

The plasma discharge occupies the region shown on figure 1. The region considered 
herein for analysis is between the potential maximum and the anode sheath. (Except for 
the absence of a magnetic field, this approximates the cathode-anode region in the ion 
thruster. ) In this region the potential is decreasing monotonically towards the anode. 
Potential defined in this manner is thus negative. A low pressure plasma is assumed, 
so that plasma ions have an effective collisional mean free path greatly in excess of the 
cathode-anode spacing. Ion motion in the plasma is determined solely by the electric 
fields developed from slight deviations from charge neutrality. These fields are 
assumed to be so small that they provide a negligible perturbation on the primary elec- 
tron energy. (Conditions for plasma shielding by space charge sheaths have been given 
by Friedman and Levi (ref. 12). ) The anode acts as a charged particle sink, so that no 
reflection of charged particles back into the discharge plasma occurs. 

In this study only the steady-state properties of such discharges are considered. In 
particular, only the plasma approximation in which the Debye length is much less 
than the characteristic plasma length L (proportional to the electron mean free path for 
ionization) is solved. 
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THEORY 


The starting point for the calculation is the Poisson equation 



■-S- (n< - n ) 
e o 


(1) 


(All symbols defined in appendix A. All quantities are in SI units. ) Following Tonks and 
Langmuir, an expression for the ion density at z (z ^ 0) is obtained as follows: Con- 
sider all positions z' < z at which ions are created by electron bombardment. It is 
desired to obtain the ion density at z resulting from ion production at all positions z\ 
This density is 


2 

njz) = f N- (z')[v(z, z ')] -1 dz’ 
• / 0 


where N^(z') is the volume ion production rate, and 


v(z, z') = j— e|V(z') - 



( 2 ) 


(3) 


is the speed of ions produced at z* arriving at z. The electron density is assumed to 
be a superposition of Maxwellian and monoenergetic (primary) components (analogous to 
ref. 2), 


n e (z) = n m (0) exp 


eV(z) 

kT 


+ n p (0) exp(-z/X) 


(4) 


where X is the total collisional mean free path for primary electrons, T is the Max- 
wellian electron temperature, n m (0) is the Maxwellian electron density at z = 0, and 
n (0) is the primary electron density at z = 0. Also note that V(0) = 0 by definition. 

r 

The ion generation is given by 


Nj(z') = Ve(z') (5) 

where 

n g(z) - n m (0) H <p exp ~j + exp | ( 6 ) 

The parameter v m is the ionization frequency for Maxwellian electrons, £ = u^/ v m is 
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the ratio of primary electron ionization frequency to Maxwellian electron ionization fre- 
quency, and <p is the electron density ratio n p (0)/ n m (0) . For mercury |, the ioniza- 

3 

tion frequency ratio, is of the order of 10 to 10 , but in cesium, £ is of the order of 
10 to 10 2 (ref. 13). 

The introduction of the dimensionless variables 


V = - 


eV 

kT 


s 



s’ = 


z* 

L 


(V 


into the Poisson equation (1) results in 


0 = 



+ < pe 


-P s 



n* (s')ds' 
V?7(s) - 77 (s’) 


( 8 ) 


T 2 nl/2 

Here Xj-j is the Debye length le Q kT/e n m (0)J , L is a characteristic length to be de- 
fined, and p is the mean free path ratio L/x. Since L in equation (8) is an adjustable 
constant, it is convenient to define it so as to render the coefficient of the integral equal 
to unity. Thus 


1 

L 



(9) 


or, alternatively, L can be defined as the reduced mean free path for ionization by Max- 
wellian electrons 


L 



(10) 


In the form of equation (8), the Poisson equation is the complete Tonks- Langmuir 
plasma- sheath equation valid throughout the plasma and sheath regions of the discharge. 
A general closed form solution to this integro-differential equation is not attainable. 
Numerical solutions for the case (p - 0 have been given by Self (ref. 6). A general 
simplification to be used in this study is that within the main body of plasma, X^ « L 
in equation (8). That is, in the plasma discharge region of figure 1, the Debye length is 
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presumed to be much less than the characteristic length L. This leads to the plasma 
condition 


n i 


n. 


(ID 


and is equivalent to the statement that plasma potential gradient changes little over a 
Debye length. The complete equation (8) thus reduces to the integral equation 


-V 


+ cpe 


-ps 


f g(s')[r?(s) - T7(s’)r 1/2 ds' = 0 

Jo 


( 12 ) 


where g(s) is defined as n*(s)/n m (0). Equation (12) states that both Maxwellian and 
primary electrons contribute to space charge neutralization as well as ion production in 
the plasma. Solutions to the plasma equation (12) for the case cp = 0 have been given in 
references 2 and 3. Solutions for the case p = 0, g(s) = 1 have been given in refer- 
ence 2. For the problem considered herein equation (12) was solved on a digital com- 
puter in the low pressure approximation ps « 1. (J. S. Sovey and R. W. Palmer assist- 
ed in performing some of the numerical calculations. ) In this general low pressure case, 
the primary and Maxwellian electrons contribute to both the ionization and space charge. 
In a later section, the limiting case in which the primaries do not contribute to space 
charge are considered. 


GENERAL LOW PRESSURE SOLUTION 

Equation (12) can be solved in the low pressure approximation ps « 1 to yield the 
plasma potential distribution t] - rj(s) (appendix B). A series solution of the form 

s = E a M x/i (13) 

P=1 


was used, where x = Vrj. Recurrence relations for the coefficients are 


l 2p+l 


2/U-l 


2(p-l) M-l 


fp^ V (q + 1 ) a q + l a 2p-q + P * E E (q+1) Vl,A + l a 2p-q-l 

1 q=0 q=0 v= 0 


P-1 


j’ 1 ' E; ‘ f E C P^ a 2p+l-2i/ - 2p(p E B P-1, ^ Ma 2 

(Zp - ip . Z l/=1 2^=0 


P 


ir(2ii + 1 ){&<p + 1) 

( 14 ) 
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where 


M = 1,2, . . . 


with 


a l " 


7< 1 + *> 

&<P + 1) 


(15) 


and 


l 2/i+2 


2p-l 


2m 


P E (q+1)A MVl a 2M-q + P ^ E (q+1) Vl a 2M+l-q 

q=0 q=0 


M 


-pcpA^Qp + l)a 2jU+1 - V d m v *2p+2-2v 

v=l 


2(ju + !)(£<? + 1) 


(16) 


where 


with 


M - 1 , 2, . . . 


pep 


(^ a l - a l) 


2 &<P + 1) 


(17) 


In the relations (14) and (16), we have 


B 


_ (-D 


M 


pv 


>2m 


L _ 1 (2 M): 

M 2 2 "M 2 
(-1)"(2m + 1) 


(2ju + 1 - 2u)(2p + 


H ‘ 

1 - t>)!z/_ 


>- 


C (2/i + 1 - 2 1/) 

v\ 


* 1 - "> 


(18) 
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Forty terms of the series (13) were computed and found to converge rapidly for rj s 1. 2 
independent of the values of 4, <P , and p used. Several solutions are shown graphi- 
cally in figure 2 for various values of the parameters. An interesting feature is noted 
from the curves of figure 2. The Tonks- Langmuir (T-L) boundary condition (which de- 
fines the limit of validity of the plasma equation) 

— = 0 (19) 

dr] 

is not satisfied for any <p 4- 0. In fact, for very small <p (<p < 10 ) and i<p < 10 , 

the solutions approach a limiting curve which closely approximates the <p = 0 solution 
for rj ~ jj. The lack of agreement for r? > ^ implies that the inclusion of a monoener- 
getic electron space charge density precludes locating the sheath boundary by solving 
the Tonks -Langmuir plasma equation. The reason for this behavior can be explained 
as follows. Transforming equation (12) by means of the Schlomilch inversion (appen- 
dix B) yields 


dr 7 _ {jcpjf ps + e _7? ) 

ds l + <p - 2Vn e'^DiVv) - p<pAv) 

where D(v'f?) is the Dawson function. For p - 0.854, where 

1 - 2Vn e~^D(y/rj) = 0 

equation ( 20 ) reduces to 

dr? = ffVr? (g <?e~ ps + e -7? ) 
ds <p[ 1 - pS(v)] 

Because the denominator of equation (22) does not vanish at r? = 0. 854 for nonvanishing 
cp , the Tonks -Langmuir boundary condition (19) cannot be satisfied. 

This is true even if the denominator of equation (22) vanishes for some r? > 0. 854 
satisfying >(r?) = (1/p), since then the factor 1 - 2 y/rj e~ ri D(s/rj) would not vanish. Thus, 
as r? — 0.854 in the limit \<p — 0, the solutions depend on the values of £, <p , and p 
and on the nature of the approximations used to evaluate J{r\). 

In general, the solutions were nearly independent of p. Thus, in the plasma, the 
potential distribution is essentially dependent only on the relative number of primaries 
entering the discharge and their ionization rate, while depending only slightly on pri- 
mary electron scattering phenomena. 


(21) 


( 22 ) 
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SPECIAL CASES - NO PRIMARY ELECTRON CONTRIBUTION TO SPACE CHARGE 
Both Primaries and Maxwellians Produce Ionization, <p « 1 

As shown in appendix B, equation (12) can be transformed by means of the 
Schlomilch inversion (ref. 14) into a form amenable to computation. In cases for which 
the primary contribution to space charge neutralization is negligible (i.e. , <p « 1), the 
Schlomilch transformation yields the differential equation 

d h _ Wrj (gffe~ ps + e~^) @3) 

ds 1 - 2Vr\ e -TJ D (\Zrj) 

The condition (19) is satisfied when the denominator of (23) vanishes. This leads to the 
usual Tonks- Langmuir value of r\ as determined by Harrison and Thompson (ref. 4). 

We thus note that the absence of primary electron space charge allows the Tonks - 
Langmuir sheath boundary condition to be satisfied independent of the nature of ion gen- 
eration. Numerical solutions to (23) are plotted in figure 3. They are functions of the 
product t,(p and, to a lesser extent, of p. In fact, at a value of ~ 6, the solutions 
were found to be independent of p over two orders of magnitude variation. 

These approximate solutions can be compared to actual solutions to the complete 
plasma equation (12) in the limits, 


Z<P ~ 1 

(p « 1 

Such a comparison is given in figure 4 in which solutions to equation (23) for | = 1 

and 10 are compared with the appropriate solutions of equation (12) in which <p 0. 01 
and 0. 001. It is apparent that a solution to equation (23) is a better approximation to the 
actual solution if, in the latter case, (p - 0.001. Again, the approximate potential pro- 
file diverges from the actual profile as the Tonks- Langmuir boundary value of r\ is 
approached. 


Only Primaries Produce Ionization, <p « 1, ij <P » 1 

Another physically interesting situation results when primaries are assumed to be 
the only contributors to ionization. The Maxwellian plasma electrons account for the 
space charge neutralization of the ions. Letting 

7] = /3a (24) 
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where 


a = -- 


0 =- 


!!« 

kT 


equation (23) becomes (using £<pe " pS » e- 71 ) 


(25) 


da 

ds 



(26) 


Now s = z/L, where L = u (2kT/M)^ // ^. Introducing a new length 
* 1/2 m 

L = v (2eV /M) 7 , a new dimensionless position s* is defined such that 
p c 

s = (L /L)s . Noting that p = (L/L )p and remembering that £ = v / , equa- 
tion (26) becomes 


da _ n Va <^e p s 

ds 1 - 2 V/3a e" /3<T D(VJa) 


(27) 


with solution 


* 


s 



(28) 


Again the Tonks- Langmuir sheath boundary condition (19) is satisfied at /3a = 0. 854, 
and the solutions depend only slightly on p . Since equation (28) describes the 
potential profile of an approximate model of the plasma discharge in which £,<p » 1 
and <p « 1, the validity of the model can be estimated by comparing these profiles with 
profiles obtained from solving equation (12). Results of such a comparison are shown in 

* _ O 

figure 5. Here solutions to equation (27) for values cp = 0. 1 and 0. 01 and L /L = 10 
are compared with actual potential profiles in which £cp =10, 100, and 1,000 and 

_f _o _3 

<p = 10 ,10 , and 10 . Agreement between the approximate and exact curves was 

_3 

best for cp - 10 , since here the condition <p « 1 is best satisfied. A solution of 

equation (23) for £(p = 10 and p = 1 is also given in figure 5. Comparison of this 
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solution with the curve solution to equation (27), shows that solutions to the two approxi- 
mate models of the discharge differ for | (p = 10 by about 6. 4 percent at the Tonks- 
Langmuir boundary, 77 = 0. 854. For values of 77 less than the Tonks- Langmuir 
boundary value, the difference is less. 


Condition for Existence of Solutions, <p « 1, 4 <P » 1 
An interesting feature of equation (28) is that real, stable solutions require that 



This situation leads to conditions on the plasma parameters. From (26) one obtains the 
condition 


/ 3 I /2 > Ml e - 0 a D (30) 

TUP 

which holds for all /3cr . In particular, condition (30) must hold for the largest value of 
2/ 7r e”^ CT D^3cr) = 1/2. 92. Thus, when substituting for /3 and p* , equation (30) be- 
comes 

/M 1 / 2 >/m\ 1 / 2 ^l E ^_ (31) 

\ T / \M / X 2.92 (p 

where T is in electron volts. The total primary electron mean free path can be written 
as 


- = — + — (32) 

X X. p X* 

where (x')~* represents the sum of the reciprocals of mean free paths for all primary 
collisions other than ionizing. Substitution of (29) into (28) yields 


a ip < 


2.92 tp 




(33) 


Relation (33) is essentially two conditions. First, it represents a condition on the ioni- 
zation mean free path when primary electron collisions other than ionizing collisions 
occur. Second, physically useful solutions require that 
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(34) 


<P> 


1 /mT \ 1//2 
2.92 \MV C / 


or, in terms of ion atomic mass A, 


(p > 


i / t y /2 

125.3 \V C A/ 


(35) 


Condition (35) holds even if no primary electron collisions other than ionizing collisions 
occur. Thus, for example, if a mercury plasma is used, relation (35) results in 


(p > 5. 65X10" 4 



(36) 


Relation (33) is thus a lower bound on the allowed value of the primary electron fraction 
in order that the approximate model of the plasma discharge described by equation (27) 
be valid. Physically, this means that strict charge neutrality in the plasma places a 
lower limit on the primary electron fraction required for ion production. An upper 
bound on (p can be estimated from figure 5. This comparison between the approximate 
solution given by equation (28) and the actual potential profiles yields, when combined 
with relation (35), 


— — < cp < 0. 01 (37) 

125.3 \V C A/ 

The lower bound in relation (37) is the condition for the existence of real solutions 
to this approximate model. The upper bound represents the limit on cp below which the 
approximate solution and the exact solution are in close agreement. 


ION CURRENT DENSITY 

The ion current density in the plasma can be calculated from the one -dimensional 
continuity equation 


— = ™*e (38) 

dz e 


where v is the ionization rate and n g is given by equation (6). 


Equation (38) states 
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that the ion current flux is balanced by the rate of ion production. Introducing the 
dimensionless variable 


s =• 


L 


where 


L - 



yields for (38) 



Defining the dimensionless current density 

I _ 1 / M X 1 / 2 J + 

+ e \2kT / n m (0) 

and integrating (39), one obtains 

n *(y') ds « 

-S 5“_ dr?’ 

n m (°) dT ?’ 



(39) 


(40) 


(41) 


For those cases in which the primary electron space charge contribution can be 
neglected, <p « 1, equation (41) yields the same dependence of ion flux on plasma poten- 
tial, independent of ion generation mechanisms. 

I + = - e'^DiVy) (42) 

71 

At the Tonks- Langmuir boundary 17 = 0. 854, equation (39) results in I + = 0. 344. This 
lack of dependence of I on the mode of ion generation is due to the nature of the plasma 
approximation. The electron density sets the ion density so as to preserve charge neu- 
trality in the plasma. Thus, as pointed out by Auer (ref. 7), the electron distribution 
in the plasma sets the ion current to the sheath. To demonstrate the role of the electron 
distribution, the ion flux resulting from inclusion of primary electrons in the space 
charge electron density can be calculated numerically from the relation 
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( 43 ) 


where 


I + = — V 1/2 + - e _?? D(VW) - ^ C S(v')dv' 

TT Tt n 



(44) 


Using a series approximation for J(r}), the resulting current density profiles are shown 
in figure 6. The major difference between equations (42) and (43) was found to be the 
first term in (43). Variations in the parameters £ and p had only a slight effect on 
the current density. Inclusion of the first and third terms of equation (43), which are 
due to nonvanishing cp , increased the current level above that which would be ob- 
tained if cp - 0. However, as in the case of the <p = 0 curve, a maximum occurs 
at some value of r,. In the case of the (f = 0 solution, this maximum occurs at the 
Tonks- Langmuir boundary value of r\ - 0.854. If we postulate that the maximum ion 
current is delivered to the sheath boundary, curves of I + (r]) would yield the values of r\ 
at such boundaries, and thus serve to spatially locate the plasma-sheath boundary. This 
postulate is justified physically since the solution for the ion current density including 
the sheath would have no maximum. These maxima in I (rj) were found to shift to values 

_Q + O 

of r/ > 0.854 with increasing cp . Thus at cp = 10 and E, = 10 , the maximum oc- 

2 

curred at 77 « 0. 854, but at cp - 0. 1 and £ = 10 the maximum was at rj « 1. 0. 


BOHM SHEATH CRITERION 

The stable sheath criterion as originally derived by Bohm (ref. 15) is a statement 
on the lower limit of the ion energy delivered from the plasma to a sheath boundary. 

m(v? ) ^ kT (45) 

where T is the electron temperature in the plasma. The criterion assumes that the 
electron space charge density is completely thermalized at temperature T. Following 
Harrison and Thompson (ref. 4), the Poisson equation in the sheath can be written 

f(v)dv - n e (46) 


=O kT 3 2 i 7 


dip 


f 


where 


14 


and 


n. 



n e = n p (0)e -ps + n m (0)e -1? 

Assuming no ion generation or ion momentum transfer collisions in the sheath region, 
ion current continuity requires that 

vf (v)dv = v b f (v b )dv b (47 ) 

1/2 

where the subscript b denotes sheath boundary. Writing v = v g 77 ' , where 
v c = (2kT/M)^ 2 , the ion density in the sheath region is 

o 

„ . ,48) 

J ^ 

Let 


F(rj) = — J v b T) -1 / 2 f(v b )dv b - n p (0)e" ps - n m (0)e v (49) 

v s 


Expanding F(?)) about r\ - 77 b yields for (46), to first order, 


e o kT 


e 


2 



(50) 


where F(rj. ) = (n- - n ) = 0 in the plasma approximation. Nonoscillatory solutions to 

o 1 e b 

equation (50) require that 


or 
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where n = n (0)e~ ps is the primary electron density and n is the total electron or 

r r ^ 

plasma density. Assuming that at the sheath edge 



one obtains the stable sheath criterion 


M<v£> a kT 





(52) 


If no primary electrons reach the sheath (n^ = 0), then one obtains the usual Bohm cri- 
terion (45). 

It is also of interest to note that in the limit p — 0, the criterion (49) becomes 


M<v£> > kT (53) 

V n m/ b 

where n m = n m (0) e_?7 and n ^ is the plasma density. This result is identical to an 
expression used recently by Masek (ref. 16) to describe conditions at the anode sheath 
in the discharge chamber of bombardment thrusters. 

Relation (52) yields a condition that replaces the Tonks- Langmuir boundary condi- 
tion given by equation (19). If primary electrons are presumed to arrive at the sheath 
boundary, then in order that the Bohm sheath criterion, given by inequality (45), be 
satisfied, one has the condition (from eq. (52)) 


dri p 

at the plasma- sheath boundary. Condition (53) states that the potential gradient at the 
sheath boundary depends on primary electron collision phenomena. It is a less restric- 
tive condition than the Tonks -Langmuir condition (eq. (19)) in the sense that the precise 
location of the sheath boundary is not defined. In the limit of very large p, relation (54) 
approaches the Tonks- Langmuir condition. 


THE MEAN ION ENERGY 

The mean ion energy distribution resulting from a monoenergetic electron density 
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distribution in the plasma can be calculated from the ion distribution function. For the 
case <p = 0, the mean square ion speed is given by (ref. 4) 


<v 2>,kT(e^-l) (55) 

1 M 

This relation is, in general, valid provided that the spatial variation of electron density 
is given by n m (0) e ”^- For the case in which np, the primary density, does not vanish, 
the plasma equation (12) can be written 

n m (0)e_77 + n D (°) e ' PS ' (— ) 1/2 L r gfa’)fa " V')~ 1/2 dr?' - 0 (56) 

m p V 2kT / L dri' 


where 


g(r)') = PS + e V ] 


The ion distribution function is given by 


f(v)dv - iJJL) g (77 ' ) (17 - t?T 1//2 — dr}' 


A Schlomilch transformation of equation (56) allows equation (57) to be written as 


f(v)dv = 


(0) 


1 + <P _ 2e~ 77 'D(v / V) _ p(f) >(77 ’) 
V (77 - V77 - 77' V?7 - 7 / 


Noting that 


f f(v)dv = 


one obtains 


n i = n m (0) (1 + <P) ~ 


n m (0) 


/ 


e~ 7j ’p(\/^ 7 ) 


77 - 7 }' 


/ 


Av'W 


V ~ V 


The mean square ion speed is given by 


f 2 ) = — y*v 2 f(v)dv 
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where 


V 2 . (W) fc - 

V M / 


) 


( 61 ) 


Substituting (58) and (61) into (60) yields the mean square speed as a function of rj, 

> 7 ] 


2 n m (0) 

<v?>=- 


/ 2kT \ 


7rn. \M / 


/ 


77 e'^m^dT t + L v<p - pep 

./ ; 2 

YT) - T) 


f 


Vv - if S(rj')dri' 


(62) 


In appendix C, it is shown that the last integrations in equations (59) and (62) can be re- 
duced, so that equation (59) becomes 


n i = n m (0) (1 + <P) - 


n rv> (0) 

m v ' 


i 


T? e~ 7?, D(V^) dr ? ; + ffy(1 - e -p S) 
Vrj - if 


(63) 


and equation (62) becomes 


2 n m ( °) 

<vf> =^— 


2k 
M ) 


r /*t? _ 

/*- 


e~ 77, D(v/ri T )d7? , + n ^ 

2 

V - rr */ 0 




/ 


r ps ' do- 


using equation (63) to eliminate the first integral in equation (62) and noting that 
n i = n m (0)(^e' ps + e _T? ) yields, finally, 


(64) 


<v 2 > = ^ 
1 M 


1 + <p 


W 


V -os' 
e ps dr]' 


<pe' ps + e”^ 


(1 - 



(65) 


For purposes of calculation, this result can be simplified in the low-pressure approxi- 
mation ps « 1 by using the substitution e~ ps « 1, so that equation (65) becomes 



1 + <p(l + 7 ?) e V _ 

epe p + 1 


1 


( 66 ) 


In the limit <p - 0, equations (65) and (66) reduce to equation (55). Equation (66) is 
plotted in figure 7 for several values of <p, including <p = 0. It is seen that the mean 
square ion speed (proportional to mean ion energy) decreases with increasing (p . In the 
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limit of very large <p , equation (66) yields the result that 

-M. <vf) « 1 J (67) 

kT 

1 2 

Thus, for large <p the mean ion energy (Ej) =g M(v. ) is given by 

< E i> “-f eV plasma < 68 > 

where Vpi asma > the plasma potential, is negative. Physically, relation (68) means 
that the average ion energy at any point in the plasma is given by the average potential 
between z = 0 and that point. A large primary electron fraction means that ion produc- 
tion is more uniformly distributed throughout the plasma than for the <p = 0 case. Be- 
cause ions are produced with a negligibly small initial kinetic energy, the mean ion 
energy at any point in the discharge would be less than that for the cp - 0 case, in 
which ion production is greatest near the rj = 0 point (where n g = n m (0)), and ions gain 
kinetic as they fall through potential 77 . 

Also evident from figure 7 is that the Bohm criterion (relation (45)) is satisfied at 
values of 77 increasing with increasing <p. For very large values of <p , relation (45) 
requires that, at a plasma-sheath boundary, tj s 1. Thus comparing such calculations 
of the mean ion energy with calculations of the maximum ion current, I , could serve 
to yield the plasma- sheath boundary value of 77 for values of <p > 0. 


PLASMA DENSITY PROFILES 

Often it is the plasma electron density that is experimentally determined from diag- 
nostic techniques such as plasma oscillation measurements (refs. 17 and 18). In the 
present problem, density profiles can be obtained from equation (4) and the curves of 
figure 2. Figure 8 gives the plasma density (normalized to the Maxwellian electron 
density at zero potential) as a function of 77 and s for two values of <p . The value of 

s is normalized by the K-factor, where K = 100 for <p =0.01 and K = 10 for 
-3 

(p = 10 . The curves show an expected decline in plasma density with increasing 77 

and s. 


CONCLUSIONS 

The results of including primary electrons in the plasma electron density in a one- 
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dimensional, low pressure plasma are as follows: 

1. The potential distribution was found to depend on ion generation mechanisms in 
the plasma, the ratio of primary electron to Maxwellian electron densities at the poten- 
tial maximum, and the ratio of primary to Maxwellian electron ionization frequencies. 

2. With the inclusion of primary electron space-charge density, the Tonks- 
Langmuir boundary condition was not satisfied for any value of the primary electron to 
Maxwellian electron ratio. An examination of the Bohm sheath criterion led to the re- 
placement of the Tonks- Langmuir boundary condition by one relating the plasma poten- 
tial gradient to the primary electron collision mean free path. 

3. The effect of primary electron scattering was much less pronounced and often 
negligible over a two order of magnitude variation in the ratio of characteristic plasma 
length (proportional to ionization mean free path) to beam electron total collision mean 
free path. 

4. Two approximate models of the plasma were considered to correspond to the 
condition of negligible primary electron fraction. The plasma potential solutions ob- 
tained for these models coincided with the exact solution except at values of potential 
near the Tonks- Langmuir boundary value. For the case in which the product of the two 
ratios of primary electron density to Maxwellian electron density and primary to Max- 
wellian electron ionization frequencies was very large, upper and lower bounds in the 
beam fraction were obtained. The lower bound is essentially an existence condition, 
whereas the upper bound was obtained by matching the approximate solution to the exact 
solution. 

5. The ion current density distribution in the plasma was also affected by the pres- 
ence of primary electron space charge. The deviation from the case of no primary elec 
trons depended primarily on the primary fraction. In the absence of a primary electron 
space charge contribution the variation of ion current density with plasma potential was 
independent of ion generation mechanisms. The variations of the ion current density 
with plasma potential exhibited maxima for different values of primary electron fraction 
By analogy to the zero beam solution, these maxima could serve to define a plasma- 
sheath boundary, assuming that the maximum ion current is delivered to the boundary. 

6. The mean ion energy in the plasma was found to be modified by the presence of 
primary electrons such that the mean ion energy at a given value of potential was less 
than the case of zero primary electrons. This result could be understood from the fact 
that the presence of primary electrons yields a more uniform ion production throughout 
the plasma. 

It is of interest to conjecture on the applicability of the preceding analysis to the 
discharge plasma of bombardment thrusters. It is apparent that primary electrons in 
the discharge play an important role in determining plasma behavior. For example, 
energy losses due to ions dissipating their energy in wall collisions could be minimized 
in the presence of a high primary fraction, because the mean ion energy would be lower 
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(see (6)) and ion production would be more uniformly distributed throughout the dis- 
charge. It is obvious, however, that the analysis is a gross simplification of the actual 
discharge environment of such thrusters. If possible, a more refined analysis should 
include the effect of a magnetic field on the electron distribution. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, October 9, 1970, 

120-26. 
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APPENDIX A 


SYMBOLS 


A 

A jU> B /I V’ C 11V’ D HV 

a 

D(x) 


ion mass, atomic mass units 
coefficients defined by relations (18) 
coefficient in series expansion, equation (13) 

,_2 

the Dawson function, 


/ 


e x dx 


d 

E i 

e 

g(s) 

K 

S(v) 

J + 

K 

k 

L 


modified Debye length (eq. (21)) 

ion kinetic energy 

electron charge 

the ratio n*(s)/n m (0) 

normalized ion current density (eq. (36)) 

integral defined by eq. (37) 

ion current density 

normalization factor in fig. 8 

Boltzmann constant 


plasma characteristic length, u 



1/2 


L 

M 

m 


plasma characteristic length, 


P 




ion mass 
electron mass 
volume ion production rate 
electron density 


* 

n g defined by eq. (6) 

n A ion density 

n m (0) Maxwellian electron density at z = 0 

np primary electron density 
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iipj total electron or plasma density 
n(0) primary electron density at z = 0 

ir 

s,s' dimensionless position variables defined in eq. (7) 
s* dimensionless position variable defined in eq. (28) 

T plasma electron temperature 

V plasma potential 

V c cathode sheath drop, V 

2 

( v. ) mean square ion speed 

z, z' position variables 

/3 eV c /kT 

e Q free space permittivity 

T] dimensionless potential variable defined in eq. (7) 

rj Q value of rj at Tonks- Langmuir boundary 

6 angle parameter used in Schlomilch transformation of equation (12) 

X total primary electron collisional mean free path 

r , 2 i 1 / 2 

X D Debye length, |e 0 kT/e n m (0)J 

/V m Maxwell electron ionization mean free path 

A. primary electron ionization mean free path in eq. (21) 

M, v summation indices 

i^m p Maxwellian, primary electron ionization rates, sec”^ 

I the ratio ^p/^ m 

p,p the ratio h/X and L /X 

u dimensionless variable, -V/V c 

(p the ratio n p (0)/n m (0) 
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APPENDIX B 


SERIES SOLUTION TO THE PLASMA EQUATION (12) 


The equation to be solved is 


•/ 


|_ ^e~ ps * + e^ ds , _ e ~rj _ ( p i 

Vrj(s) - T7'(s') 


-ps 


(Bl) 


Let 


Then 


2 , 2 . 2 n ,2 

T) - x , * = x sin 0 = x’ 


V?7 - ?]* - x cos 0 = ^2- 
d0 


and equation (Bl) transforms to 

nr/2 


/ n/t , 

V 


5¥ , e -ps' + e -x' 2 ) dsl de . e -x 2 _ ?e -ps(x) 
1 dx’ 


(B2) 


The Schlomilch inversion formulas (ref. 14) can now be applied. 

2 


f(x) = e x + cpe 


■ps 


^M = -2xe" x -pcp^-e-PS 
dx dx 


h(x’) =- (| <pe~ ps ' + e _x ' 2 ) — 
o V / dx’ 


(B3) 

(B4) 

(B5) 


We note that 


2 n ^ 2 
f (x) - — J Q h(x’)d0 


(B6) 


where 
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r n/2 

h(x) = f(0) + x Jq f'(x sin 0)d0 


(B7) 


where the prime on f denotes differentiation with respect to its argument. Substituting 
(B3) to (B5) into (B7) yields 


— cpe~ ps + e -x2 ) — = 1 + cp - 2x 


dx 


/ 


ff /2 o . 2 

e“ x s * n 0 x sin d d0 


-pcpx 


L 


n/2 


e-ps’dslde 

dx' 


9 o 

Using the identity sin = 1 - cos allows equation (B8) to be written as 


p<px 


L 


ff//2 e" ps ’ — d0 = 1 + (p - 2xe x D(x) 
dx' 


_lf^ e - ps + e- x2 )^ 

2 ' ' dx 


(B8) 


(B9) 


Equation (B9) is thus the transformed plasma equation (12) of the text. 
To equation (B9) the series 


s = Z a n xl1 

n=l 


(BIO) 


is applied. The factor e ps requires an approximation. For low pressure plasmas 
ps « 1, so that e” ps « 1 - ps. Then 


,-ps ds 





dx 


‘-*Z 

k=l 


a k x 


S (q + 1)a q+l xq 

Lq=o 


(B11) 


and 


/-»/2 ds , a . r n 

f e P ^7 de “ Z <1 * »Vl* q / 

y 0 q=0 y 0 


sin q 0 dd 


w w jr*jr /O 

- p E Z (q + l)a kVi xkxq f sin +q0 de (Bl2) 

~ n i. i 0 


q=0 k=l 
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In order to evaluate the integrals on the right side of equation (B12), the sums need to 
be separated into even and odd powers of x. 


S 1 = S (2/x + 1)a 2ju+l x2M / 7 sin2M0 dd 

M=0 0 


M=0 u 


d0 


From tables (ref. 19) 


rn/2 

Jo 


sin 2 ^0 d0 


1_ /2 ii\ n 
2u V u / 2 


2 2ju Vpi 


and 




jt/2 


M 


sin 2M+1 0 d0 = 


J_ (.dm'P Hi! /2 m + 1\ 

2 ii / j 2[i + 1 - 2v \ v / 


v=0 


so that 


oo 

-■D 


S, = > 


>2 /i 




ju=0 


OO 

■E 


(- 1 ) 


M 


+ 2 7 .Err + 1)a 2M+2 


,2/i 


li=0 


Ii 


E 

i/=0 


ME 


(2M + 1)! 


(2m + 1 - 2v)(2li + 1 - v)\v\ 


2fi+l 


For the second sum, 


S 2 - P 


OO OO 

EE 


<q + l>a k a q+1 x k+ <l / ^ sin^O de 




q=0 k=l 

One must account separately for k + q even and odd. 


(B13) 


(B14) 


(B15) 


(B16) 


(B17) 
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co 2fl -1 

S 2 = P (q + l)a, 

ju=l q=0 


q+1 


“^-q^ X 


jt/2 


sin 2p 0 d0 


°° 2p /2 
+ p ^ ^ (q + l)a q+1 a 2(J+1 _ q x 2 ' 1+1 sm 2l ‘ +1 e d8 (B18) 

p= 0 q=0 

Substituting equation (B14) and (B15) into (B18) yields the series form for S 2 - Thus, 
the integral (eq. (B12)) can be written 


/ jt/ 2 00 °° ju 

e ps ^s_ d0 ^ (2 m + + 2 2 2 B M> + 1)a 2p+2 x 

^ ^ p=0 p=0 v=0 


,2 P+1 


00 2p-l 

~1 P Yj E (q+ 1)A M a q + l a 2p-q x: 

z p=l q=0 


2p 


00 2p p 

P Z 2 E (q + 1)B Mi' a q+l a 2n+l-q X 

p = 0 q=0 i^=0 


2 p+ 1 


(B19) 


where 


A - 1 ( 2 P) ! 

M 2 2m (p!) 2 


(B20) 


and 


B 


liv 


_ (-D M 

>2p 


(-i)^(2p + l): 


(2p + 1 - 2iv)(2p + 1 - 


(B21) 


The series expansions of other terms in equation (B9) are 


00 00 

e -ps * £ (2p + l)a 2M+1 x 2M + 2 £ (p + l)a 2jU+2 x 

OX , , r» . . _n 


2 p+i 


p = 0 


p=0 

2p 


to 2p-l 00 2p 

" P E E (q+ 1)a q + l a 2p-q x ^ " P E E (q+ 1)a q + l a 2p + l-q x 

p=l q=0 p=0 q=0 


2p+l 


(B22) 
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(B23) 


2xe~ x D(x) 


CO 



M= 1 


(-1) m 2 ^x 2m 
(2ju - 1)!!. 





^ = E E C M ^ a 2 M + l - 2 ^' 

^ /u=0 i/=0 


2 m 



V D a x 2/i+1 

2_J u i±v' a 2M+2-2i/ x 

i^-0 


(B24) 


where 


C l(1 <2ju + 1 - 2v) (B25) 

M v v'. 

and 

D = 2 il?E (ju + 1 - v ) (B26) 

vl 

Now, we substitute equations (B19) and (B22) to (B24) into (B9) to obtain the series rep- 
resentation of the plasma equation. 


§ M E V 2/J + 1)a 2M + l x2M+1 + 2 ^El B M> + 1 > a 2 # x+2* 2#i+2 

M=0 M=0 t/=0 


00 2/J.-1 co 2m m 

1 p * E E (q + 1)A M a q + l a 2M-q x2,i+1 ' ^ E E £ (q + ‘Vqd^^l-q* 2 " 

^ = 1 q-0 m =0 q-0 n=0 


+2 


QQ 


(2m - 1)!! 2 


TT + i ^ E (2M + 1)a 2M+l x2M + *(v E (/i+ 1)a 2M + 2 x2M+1 




/I-O 


4=0 


00 2/i-l 2/i 

- j Pi(fi E E (q + 1) Vl a 2M-q x2M - f E E (q + 1)a q + l a 2M + l-q x2M+1 


°° 2/i 


M=1 q=0 

OO /i 


/i=0 q=0 


2 EE C Mi' a 2M+l-2i' X M + 2 E E D Mi' a 2M+2-2^ x2A1+1 " 1 ~ ^ " 0 (B27) 

/i=0 i/=0 /i-0 p'-O 
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The expansion (eq. (B27)) is then used to obtain the recurrence relations in the coeffi- 
cients a n of the series (eq. (BIO)). 

For /i = 0, one obtains 


2^ a 1 + 2 C 00 a l = 1 + <P 

and since Cqq = 1 

_!»♦*> 

1 ~~ iv 7 T 

Also, from the coefficient of x\ 

W- A Qai + H <P&2 ~ \ + l D 00 a 2 = 0 


(B28) 


(B29) 


(B30) 


so that 


p<p (|a^ - a^) 
2 &<p + 1) 


(B31) 


For ju > 0, one obtains 


1 2m+2 


2(M + 1 )(k<P + 1) 


2 / 1-1 


2ju 


p * S (q+ 1)A M a q+l a 2M-q + P ^ E (q+ 1} Vl a 2M+l-q 

q=0 q=0 


M 


- p V A (j (2m + l)a 2t , +1 - £ D M1/ a 2(J+2 . 2j/ 

V=1 

For the coefficients of odd powers of x, and from even powers of x 


(B32) 


l 2ii+l 


2/i-l 2(/i-l) /i - 1 


q=0 


q=0 1^=0 


(-1) M 2 P 7 T 


fcLi 1 


TT^'lE C M, a 2M + l-2,- 2p<p Z MB M-l,^ a : 


(2p - 1)!! 2 ^ 


^=0 


2m 


7t(2m + !)(£<? + 1) 


(B33) 
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Using equations (B29) and (B31) to (B33) the series (BIO) was computed for 40 terms. 
The parameter | was varied from 0. 1 to 100, <p was varied from 0. 00016 to 0. 5, and 
p was varied from 0. 01 to 1. 0. The convergence of the series was quite rapid for 
fi s 10 and 77 s 1.21. 
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APPENDIX C 


REDUCTION OF CERTAIN INTEGRALS 

In the calculation of mean ion energy in the plasma, two integrals occur that require 
special treatment. Since the procedure for reducing these integrals is identical, a de- 
tailed treatment will be given only for one of the integrals. 

Consider the integral occurring in equation ( 56 ) 



>(T)')dr?' 

Vv - V~' 



f(77”)dT?" d n' 

Vv - t)' vV - i?” 


(Cl) 


where 

f(Tj) = e~P s 

dri 

In order to reverse the order of integration, we introduce the step function 

CO for rj M > n' 

A (17* - 77") = j 

ll for 77" < 77' 

to eliminate the 77' upper limit in equation (Cl). Thus 



A ( 77 ’ - 77 ,, )f( 77 ,, )d7]" dr)’ 
^77 - 77 ' V^T)' - 7]" 



Vv ~ 77 * 


Let 


I = 



A (7 7' - 77 T 1 )drj * 
V77 - 77' V 1 ?' “ T ?'* 


The following sketch can be used to evaluate equation (C 5 ). 


(C 2 ) 


(C 3 ) 


(C 4 ) 


(C 5 ) 
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For rj" < Tj, equation (C5) yields 


1 = 2 tan 


1 

v V* - v” 


= - 71 


For 77 ” > rj, it is apparent that 1 = 0. Thus, equation (C4) becomes 



Siri'ldr]' 

Vv ~ r? 


-it 



f(T7 M )dr} M 


From equation (C2), we can write 


f(ij) (e' pS ) 

P dr? 


(C6) 


(C7) 


so that 



v SW)dr£ 

Vv - rp 


= — (1 - e" ps ) 
P 


(C8) 


The second integral of interest is, from equation (60), 



Jiv'W 


■ff&r 


f(r?")d77” d rj' 


(C9) 


Again, using the step function (C3) to replace the r ?' upper limit by °° allows reversal 
of the order of integration. In this case 



(CIO) 


so that for r/” < t] 


I = 0?-TJ")- (C11) 

2 

and, again, 1 = 0 for q" > rj. In this way the integral (eq. (C9)) becomes 
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J ( Vv - v' Siv'W =-Jq bn - 

2 


(C12) 


Integration by parts of the right side of equation (C12) yields 






(C13) 
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Normalized plasma potential, 



Dimensionless position, s 


Figure 5. - Comparison of exact solutions with approximate solutions, <p«l. (} = eV c kT;p* - Vji/£)p. 


Ionization Electron Mean free 
potential density path ratio, 
ratio, ratio, p 
€ 


100 0.1 1 

10 .1 1 

100 .01 1 



0 .2 .4 .6 .8 1.0 1.2 1.4 

Normalized plasma potential, 


Figure 6. - Normalized ion current density. 
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Figure 7. - Ion mean square speed against normalized 
plasma potential. 
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